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The optimization problems denning meta-stable or stationary equilibrium are explored. The Gibbs 
scheme is modified aiming to describe the statistical properties of a class of non-equilibrium and 
metastable states. The system is assumed to maximize the usual definition of the Entropy, subject 
to the standard constant energy and norm restrictions, plus additional constraints. The central as- 
sumption is that the existence of the considered metastable state is determined by the action of this 
£Nj | additional dynamical constraint, that blocks the evolution of the system up to its maximum En- 

tropy state, thus maintaining it in the metastable or stationary configuration. After requiring from 
the statistical description to be valid for the combination of two nearly independent subsystems, it 
| follows that the eigenvalues of the constraint operators C should have the restricted homogeneous 

form C(pi) = p\ , in terms of the eigenvalues of the density matrix pi, where q is a fixed real number. 
Therefore, the distribution of the eigenvalues of the constraint function have the Tsallis structure. 
This conclusion suggests the interpretation of the q parameter as reflecting the homogeneous de- 
pendence of the constraint determining the metastable state on the density matrix. An application 
to the plasma experiments of Huang and Driscoll is expected to be considered elsewhere in order to 
' compare the results with Tsallis scheme and the minimal Enstrophy one. 



a 

W ■ I. INTRODUCTION 

The statistical physics of systems which remains for long times in metastable and stationary states is an important 
theme of research nowadays [2||3||j|- Since the appearance of the proposals for variation of the Boltzman statics 
*^ ■ proposed by Tsallis for a wide class of non-equilibrium states, it has a been an explosion of interest in the field 

A particular point of attention in the literature is devoted to the understanding of the meaning of the special real 
parameter q which fully characterize the deviation of this new statistical description from the Boltzman one 0Q| . 

In the present work we explore the statistical properties of metastable and non-equilibrium stationary states with the 
aim of determining possible new ways of studying these properties and their possible connections with the conceptual 
elements in the Tsallis approach. Out central starting assumption is that statistical properties differing from those 
of the Gibbs thermal equilibrium states, could be generated, at least in a large class of systems, by the existence of 
a dynamically generated constraint stopping for a while their relaxation to the thermal equilibrium, described by the 
Bloch density matrix. Since the systems under consideration are in quasi-equilibrium in their metastable states, it is 
here assumed that the additional constraint is approximately conserved. Then, the procedure being proposed here, 
consists in assuring that the system on its evolution, tends to maximize the usual Entropy in terms of the density 
matrix S = —k Tr [p log p], under the standard restrictions of constant energy and norm, plus an additional specially 
constructed constraint. 

The commutation of the statistical density matrix p describing the metastable state and the Hamiltonian is assumed 
JL* ■ since the physical quantities of the metastable state are quasi time independent. Consider now the evolution of the 
system in the quasi-equilibrium state. It is clear that the probability distribution of the energy values of the Gibbs 
subsystems will not follow the one given by the Bloch density matrix. There will be a different distribution along the 
metastable evolution defined by a density matrix p. This behavior is assumed here to be produced by the existence 
of a constraint C which constant value during the movement means its conservation. That is, it will be assumed to 
commute with the Hamiltonian. Also considering that the energy spectrum is non-degenerated, the three operators 
H, p and C can be simultaneously diagonalized. This property allows to express the constraint C as a function of 
the density matrix C[p\. Taking into account that the product C[pH is also a conserved quantity, a natural way for 
imposing the effects of the constraint C in the variational Gibbs problem arises: when maximizing the Entropy the 
constant value of the trace of C[p]H will we also required. In other words, this form of the new restriction will added 
in this paper onto the standard Gibbs maximization problem defining the statistical density matrix by the conditional 
extremum point of the Entropy as a functional of the density matrix. 

As a consequence of the above definitions, the interesting result follows that the additivity of the modified conserved 
mean values, for two approximately independent subsystems, directly implies that the constraint function C[p] should 
have the Tsallis homogeneous structure C(p) = C q p q , with q being an arbitrary real number. The possibility that the 
density matrices solving the new conditional extremum problem corresponds to metastable states is furnished by the 
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fact that, by assumption, the constraint C is a weakly conserved quantity, which can be unstable because in general, 
there could no exist a rigorously conserved dynamical variable associated to it. However, it could has the chance of 
locally constraint the system for a while, to be in a metastable state. 

The equations of the extremum problem for the determination of the density matrix describing the considered 
metastable or non equilibrium stationary state are also presented. 

The work will proceed as follows. In Section 2 the proposal of the scheme is presented. Section 3 is devoted to 
argue that, if the statistical description of physical system satisfies the additivity condition, then the dependence 
of the constraint C on the probabilities have the Tsallis structure. The explicit form of the extremum equations 
determining the density matrix of the metastable state are also presented in this Section. Finally, in the summary, 
the conclusions are reviewed and some possible extensions of the work commented. 



II. STATISTICAL MECHANICS OF SOME METASTABLE AND STATIONARY STATES 

Lets us consider the quantum description of a physical systems having a dynamics fixed by a Hamiltonian operator 
H. In the Gibbs approach the properties of the systems in thermal equilibrium are described by the Bloch density 
matrix 

p = exp(-— ), (1) 

satisfying \H, p] — and which is the conditional maximum of the Entropy functional under constant mean energy. 
That is, the Bloch density matrix is the extremum of the functional 

S = -k Tr[p log(p)] + a(Tr[ P H] - E) + 
+ f3(Tr[p]-l), (2) 

in which a and (3 are auxiliary Lagrange multipliers for imposing the conservation of the energy E and the total 
probability equal to 1, and k is the Boltzman constant. 

Let us assume that during a relatively large relaxation time r, the system is not allowed to approach the thermal 
equilibrium state, because it is in a metastable or non-equilibrium stationary state. We will assume that this 
metastable state is produced by a dynamically conserved quantity C which approximately restricts the rapid evolution 
of the system to the Gibbs thermal equilibrium. The approximate satisfaction of this constraint for large times will 
be represented here by its, also approximate, commutation with the Hamiltonian. 

[H, C] = 0. 

The stationary character of the density matrix p will be assumed to also imply [H, p] = 0. As mentioned above, in 
order to simplify the discussion, let us consider that the spectrum of the Hamiltonian is not degenerate. Therefore, 
H, p and C are all diagonalized in the common basis of eigen-functions of H. Henceforth, C can be expressed as a 
certain function of the density matrix C = C[p\. 

But, due to the above definitions the quantities C[p]H is also conserved 

[C[p}H,H}=0, 

[C[p],H]=0. (3) 

Therefore, a convenient way emerges for including the effects of the constraint C in the Gibbs maximization problem 
determining an extremum of the entropy S, in the considered metastable or stationary state. The proposal is to 
include the constancy of the specially defined quantity 

Tr[C(p)H] 

Tr[C(p)} ' 1 > 



through a corresponding Lagrange multiplier in the extremum problem. Then, adding this new constraint to the 
expression of the Entropy after to be multiplied by the corresponding Lagrange multiplier, the conditional Entropy 
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functional takes the form 

S = -kTr[plog(p)} + a(Tr[pH] - E)+ 
= -kJ^Pt Iog(Pi) + a^PiCi -E)+ (3{J2pi - 1) + 

i i i 

In what follows we will assume two possibilities: 

1) The a parameter is non vanishing and becomes an effective Lagrange multiplier. This corresponds to impose 
both: the mean energy conservation constraint in common with the C dependent constraint. 

2) The a parameter is taken as vanishing and then, it is disregarded as a Lagrange multiplier. This variant 
corresponds to only fixing the conservation of the new constraint in addition to the constant probability one. 



III. TSALLIS q PARAMETER FROM ADDITIVITY OF THE DESCRIPTION 

Let us now consider the implications on the modified problem produced by requiring the additivity of the statistical 
description. That is, the condition will be imposed that the application of the scheme to a combination of two nearly 
independent systems, each of them being in the same kind of metastable state, should be equivalent to the separate 
application to each one of the systems. 

Therefore, consider two similar and weakly interacting subsystems, each of them being in a metastable state of 
the same sort. The Entropy and the constraint functions of the first system after the density matrix furnishing the 
maximum has been found, are given by 

S« = -fc^ P f ) log(^ 1) ), 

i i 

Analogously for the second system, these same quantities take the form 

sV^-k^Logip™), 

£ (2) = E^M 2 \ i = E^ 2) > 

i i 
/-./ (2)\ (2) 

c ~ E.c(pf')] ■ (7) 

The description is assumed to be also valid for the combined system. In this case the Entropy and constraints for 
the composite body are 

(id) 

£7(1,2) _ V D (1 ' 2) e (1 ' 2) 1 - V o (1,2) 

^ _ P (i d) e (i,3) ' 1 - P (id) ' 

(id) (id) 
p (l,2) _ ^(id) U[ P(id) )e (id) ^ 

^(id) U[ -P(id) > 
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where symbolizes the couple of indices defining the states of the combined system \i,j) = \i) x Then, the 

independence of the systems allows to write for the probability of the combined state \i,j) the relations 

(1.2) (1) (2) 

EW = E^E^ = 1 - ( fl ) 

1 3 

Therefore, the constant probability constraint of the individual components imply the same property for the combined 
system. 

The additivity of the Entropies follows in similar way as usual 

* (1 ' 2) = - 

(i,3) 

= -*EfW 9)l °8(Pi (1 M a) ) 

(i,3) 

= -k^p™ logfe (1) ) -k^P? I«g(pf ) 

= 5«+^ 2 ), (10) 

as well as the addition of the mean energies 

E (i,2) _ y- (1,2) (1,2) 

E(l) (2), (1) . (2V 

(i,3) 

Efl) (2) (1) . (1) (2) 



'l.„(2.,(l. _ V" „<1>,.' 2 U 2 ) 

(i)ji) ± r J 2 U 2 ) 



Now, we will also assume that the modified C dependent mean value also satisfies the statistical independence 
condition for the combined system 

where C^'^ = C(p^'^), — C{p^p) and = C{p^). Therefore, the specially defined mean value also satisfies 
the additivity properties as follows 

v r^( 1 '2) e (i,2) 

,71(1,2) _ M^j) °(i,j) e (i,j) 
- — (i,2) . 

y^ /^(l)/o(2) 
y, ^(1)^(2) (1) y, r ,{l) r ,{2) (2) 



V- ^(1)^(2) y- r (l) r (2) 

E 4 cf EiCf 

,(1) , p(2) 



4 1) +4 2) . (12) 
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However, this condition also imposes a strong restriction on the possible forms of the function C defining the modified 
mean value. In order to see this, the statistical independence condition (JTTJ) can be rewritten in the form 

C{p^pf) = C(p^)C(pf). (13) 

Considering C as expanded in powers 



n=0 



C(x) =x v Y,fn x n (14) 

and substituting this relation in (fTS)) . it follows 

OO OO CO 

(?f Vfr E /«( pfr = {^nvfr E E pPh *>? h 

n— n— m— 

which after cancelling the common factors can be rewritten as 

OO OO 

= E E m pPn p^nsnm m. (15) 

n— m— 

Next , after taking into account the completeness of the basis of powers of the variables, and assuming that the values 
of the pi variables form a continuous for all values of i, it follows 

fn(S nm ~ fm) = 0, for all m and n . (16) 

Henceforth, assuming that a particular f mo ^ 0, implies that for all m ^ m a the f rn expansion parameters are 
equal to zero. Thus, the only non- vanishing coefficient is the /„ io .one. Consequently, the allowed forms of C are 

C{ P ) = !rnJ n °+ V 

= h P q - (17) 
A. Extremum equations 

After the form of the function C has been determined, the Entropy functional as modified to impose the constraints 
through the Lagrange multipliers procedure, can be written in the more explicit form 

S = -k E Pi Iog(ft) + a(E Piti ~E) + /3(E Pi ~ !)+ 



< HiP\^ 
T,iPi 



- 1{ ^±-E q ). (18) 



where the index C in Ep has been changed to q . 
The extremum equations are 

dS „ dS „ 

— = 0, ^-=0, i = l,2...oo 

opi oa 

which after their explicit evaluation leads to the following set of coupled equations for the eigenvalues of the density 
matrix p 

9-1 

Q P 

k - a e, - (3 = -fclog(p,) + 7 ' q (e l - E q ), 

2-~ii Pi 

i 

1 = E^' 

i 

E q = (20) 
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Multiplying the first equation by pi and summing over i, a relation between the Entropy and the lagrange multipliers 
follows 

q 

= -k Y,P* 1o §(p0 + 7 E 4^1^ -E q )-l + aE + 0, 

i i Pi 

= -k^pi log^) - k + aE + f3 

i 

= S + aE -k + /3. (21) 

As usual it could be employed for constructing a generalization of the Free Energy and other Thcrmodynamical 
Potentials. However, we will delay this discussion to future extensions of the work. 

IV. SUMMARY 

A statistical description of metastable states is proposed. The main idea is that a large class of metastable and 
stationary states could be associated to the existence of quasi-conserved constraints blocking the approach of the 
system to thermal equilibrium. Assuming that this is the case, it is shown that if the statistical properties of the 
system have the additivity property, the eigenvalues of the constraint function in terms of the probabilities have the 
Tsallis structure p\. This outcome indicates the interpretation of the Tsallis parameter as furnishing the degree of 
homogeneity of the constraint as a function of the density matrix. Finally, the proposed description is defined by 
the density matrix which maximizes the usual expression of the Entropy under the conservation of the mean energy, 
normalization and the additional constraint assumed to be enforcing the metastable state. 

In future extensions of the work, it is planned to determine the predictions of the analysis for the plasma experiments 
of Huang and Driscoll Q . It is expected that the procedure could appropriately describe the experimental results of 
Ref. ■ Those results were reasonably explained by the minimization of the so-called Enstrophy function Q and also 
by a study employing the Tsallis statistics in Ref. [2J . Thus, the consideration of this task seems to be an appropriate 
step in checking the predictions of the discussion presented here. 
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